Morphology and magnetism of multifunctional nanostructured 7-Fe203 films: 

Simulation and experiments 

R.F. NeumaniQ and M. Bahiana 
Instituto de Fisica, Umversidade Federal do Rio de Janeiro, 
Caixa Postal 68528, Rio de Janeiro 21941-972, Brazil 

L.G. Paterno and M.A.G. Soler 
Universidade de Brasilia, Instituto de Fisica, Brasilia DF 70910-900, Brazil 

J. P. Sinnecker 

Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro RJ 22290-180, Brazil 

J.G. Wen 

University of Illinois at Urbana- Champaign, Frederick Seitz Materials Research Laboratory, Urbana IL 61808, USA 

P.C. Morais 

Universidade de Brasilia, Instituto de Fisica, Brasilia DF 70910-900, Brazil 
(Dated: September 1, 2011) 

This paper introduces a new approach for simulating magnetic properties of nanocomposites com- 
prising magnetic particles embedded in a non-magnetic matrix, taking into account the 3D structure 
of the system in which particles' positions correctly mimic real samples. The proposed approach 
develops a multistage simulation procedure in which the size and distribution of particles within 
the hosting matrix is firstly attained by means of the Cell Dynamic System (CDS) model. The 3D 
structure provided by the CDS step is further employed in a Monte Carlo (MC) simulation of zero- 
field-cooled/field-cooled (ZFC/FC) and magnetic hysteresis loops (Af x H curves) for the system. 
Simulations are aimed to draw a realistic picture of the as-produced ultra-thin films comprising 
maghemite nanoparticles dispersed in polyaniline. Comparison (ZFC/FC and M x H curves) be- 
tween experiments and simulations regarding the maximum of the ZFC curve (Tmax)) remanence 
(Mu/Ms) and coercivity (He) revealed the great accuracy of the multistage approach proposed 
here while providing information about the system's morphology and magnetic properties. For a 
typical sample the value we found experimentally for Tmax (54 K) was very close to the value 
provided by the simulation (53 K). For the parameters depending on the nanoparticle clustering the 
experimental values were consistently lower {Mr/Ms — 0.32 and He = 210 Oe) than the values we 
found in the simulation [Mu/Ms = 0.53 and He = 274 Oe). Indeed, the approach introduced here 
is very promising for the design of real magnetic nanocomposite samples with optimized features. 

PACS numbers: 75.70.-i, 75.75.-c, 75.75. Cd, 75.47.Lx, 02.70.Uu, 05.10. Ln 



I. INTRODUCTION 

Assemblies of nanoscaled particles or clustered nanosized systems have attracted a great deal of interest either 
due to their fundamental properties or their potential use in several technological application^^. In particular, 
composites containing nanosized magnetic particles dispersed in a hosting matrix can be envisaged as multifunctional 
materials with promising applications in emerging field^^HH, in tJ^is particular class of material system isolated 
magnetic particles are sufficiently small, so that they usually remain in a single-domain magnetic state^. Nonetheless, 
it is possible to modulate the particle-particle interaction over a broad range of magnetic particle concentration in such 
a way that inter-particle interaction starts playing a significant role on the nanocomposite end propertied'''®'. Different 
types of magnetic interactions are found in nanoparticulated systems, including the magnetic dipole coupling; the 
indirect exchange or Rudermann-Kittel-Kasuya-Yosida (RKKY) interaction mediated by conduction electrons, for 
instance when the matrix and particles are both metallic; the short-range direct exchange interaction when surfaces of 
neighbor nanoparticles are in close contact; and the tunneling exchange inte ractio n, which occurs when nanoparticles 
are considered to be only some tenths of nanometers apart from one anotheil2Eni. 

Besides the complexity of the dipolar interaction there is also an uncertainty in determining the anisotropy axis 
direction and location of the magnetic nanoparticle within the nanocomposite matrix. However, advances in theo- 
retical, experimental and simulation have contributed to unveil the magnetic behavior of one and two dimensional 
systems of interacting particleJ^^ '"^^'. Monte Carlo (MC) simulations have been used to assess the magnetic properties 
of nanocomposites comprising iron oxide-based nanosized particles embedded in non-magnetic matrices. Bujan-Niifiez 
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et alP^ studied the influence of the nanoparticle size on the blocking temperature of interacting systems (dipoledipole 
interaction) using MC simulation of zero-field-cooled curves. Hoppe et al!^ employed a similar method to study the 
influence of the dipolar interaction on the blocking temperature of a random array of superparamagnetic particles 
(maghemite) dispersed within different polymeric matrices. Lamba and AnnappornP^ reported simulation results 
including effects of disorder, short-range exchange, long-range dipolar interaction and anisotropy effects, using pa- 
rameters of a maghemite-based system. 

It becomes clear that the simulation of the magnetic response of realistic 3D systems of interacting nanoparticles is 
mandatory to understand the relationship between the nanocomposite structure and its magnetic properties. Despite 
efforts expended to address this issue not many investigations have been developed with great level of detail. For 
instance, assemblies of monodisperse particles placed at lattice positions have been used to describe experimental 
systems with much more complex morphology. Another possible approach, but rather inefficient, sorts particle sizes 
from a distribution and places them at random positions. Since particles cannot overlap, the rejection frequency for 
new positions gets higher while the sample becomes denser and, consequently, the algorithm becomes very slowly. 

In the present study we propose an alternative approach, a multistage simulation in which the 3D spatial structure 
of a nanocomposite system comprising magnetic nanoparticles dispe rsed within a non-magnetic matrix is firstly 
generated by means of the Cell Dynamic System (CDS) modePlSil. with the use of an iterative algorithm the 
model is capable of generating the main systems features, namely particle positions, sizes, and shapes. In a second 
step, the generated 3D structure is employed in the MC simulation of zero-field-cooled/field-cooled (ZFC/FC) curves 
and hysteresis loops. The simulated 3D structure and the corresponding magnetic properties are then compared 
with the experimental data recorded from real samples, namely ultra-thin films of polyaniline (PANI) encapsulating 
citrate-coated maghemite (cit-MAG) nanoparticles. Our findings successfully proved that the computational approach 
employed in the present study can be used to assess, with great accuracy, even for more complex systems, the major 
morphological and magnetic features of the real samples. A similar approach has been used in simulation of ZFC-FC 
curves of CuCo film^^^, but in this case the simulated sample was two-dimensional and the cobalt nanoparticles had 
a broad size dispersion and no local order. 

II. EXPERIMENTAL PROCEDURE 
A. Synthesis of Nanoparticle dispersions and Preparation of Nanostructured Films 

Magnetic nanocomposites were produced via the layer-by-layer (LbL) technique following the procedure previously 
reporte d ^ in which layers of anionic citrate-coated maghemite (cit-MAG) nanoparticle and cationic polyaniline 
(PANI) were alternately and electrostatically adsorbed onto a solid substrate. The cit-MAG suspension was prepared 
in a two-step chemical route as reported elsewher^23 and briefly described as follows. In the flrst step Fe(II) and Fe(III) 
ions were co-precipitated in alkaline medium to produce magnetite nanoparticle, which was oxidized to produce 
maghemite. In the second step, the as-produced maghemite nanoparticles were rendered negatively charged via 
functionalization with excess of citric acid and suspended as an aqueous-based magnetic fluid sample. The doped 
PANI solution was prepared by dissolution of a commercial undoped PANI (Mw 10,000 g mol^^, Aldrich, USA) in a 
mixture of N',N"-dimethylacetamide and hydrochloridric aqueous solution (pH 2.7). The nanocomposite production 
was realized by alternate and successive immersions of hydrophilic (lOO)-oriented Si stripes into PANI (polycation) 
and cit-MAG (anion) solutions. In between immersions substrates were rinsed in a stirred bath of HCl solution (pH 
2.7) and dried with N2 flow. Multilayered magnetic nanocomposite films of (PANI/cit-MAG)^, with £ increasing 
from 1 to 50 bilayers, were produced by repeating the cycles described above. Samples with ^ = 10, 25 and 50 were 
systematically investigated. All LbL depositions were carried out at room temperature (w 25 °C). 

B. Measurements and Instruments 

Surface morphology of the thin composite films was characterized with scanning electron microscope (SEM) (Hitachi 
S4800). Transmission electron microscopy (TEM) micrographs were recorded with a 200 kV instrument (JEOL JEM 
2100) to evaluate the average particle size and size distribution of nanoparticles in the magnetic fluid sample as well 
as their distribution within the nanocomposite films (cross-section images). Figure [ija) shows the surface morphology 
of a (PANI/cit-MAG)io film, with nano cit-MAG spheres with little size dispersion and local HCP packing. 

Values of the average particle diameter (£>tem = 7.5 ± 0.1 nm) and standard diameter deviation (tr — 0.33 ± 0.03 
nm) of the suspended magnetic particles (magnetic fluid sample) were obtained after fitting the particle size histogram 
using a log-normal distribution function, as described in the literatur^25]_ Preparation of the Si substrate-deposited 
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nanocomposite films for cross-sectional TEM observation followed the standard protocol: films were firstly glue face- 
to-face, then Si substrate was mechanically ground down to about 20 fim in thickness, following the ion milling of the 
sample to perforation. The resulting histogram can be seen in Fig. ^a,). 

Magnetic properties of the nanocomposite samples, at different applied fields and temperatures, were measured 
using a Cryogenics S600 SQUID magnetometer. The diamagnetic correction was taken into account while analyzing 
the data. For the ZFC/FC measurements the samples, as well as the measuring system magnet, were initially 
demagnetized at room temperature, ensuring the zero field condition. After demagnetization samples were cooled 
down to 2 K when the measuring field (36 Oe) was turned on. The ZFC measurements were taken while heating 
the sample at a constant rate of 1 K/s. After reaching room temperature, the sample was cooled down to 2 K 
again, without switching off the measuring magnetic field (36 Oe). FC measurements were taken while heating the 
sample at the same rate (1 K/s). Magnetization (M) versus field (H) hysteresis loops were taken at different fixed 
temperatures (2 and 300 K). The total magnetization of each nanocomposite sample was scaled to the %wt content of 
cit-MAG, the latter obtained from UV-vis spectroscopy (Shimadzu spectrophotometer, model UVPC 1600). The %wt 
contents of cit-MAG within the nanocomposite samples were obtained by measuring their optical absorption at 480 
nm wavelength. The cit-MAG contents were determined from a previously constructed calibration curve based on the 
absorbance values of cit-MAG suspensions, at different nanoparticle concentrations and at 480 nm, which corresponds 
to the typical wavelength for electronic transitions in iron oxided^. 

III. SIMULATION 

The simulation procedure comprises three main stages: generation of the spatial structure, attribution of the 
magnetic properties and simulation of magnetization curves. For each part different computational techniques were 
used, as described below. 

A. Three dimensional simulation of the samples morphology 

The (PANI/cit-MAG)f nanofilms investigated are a typical binary system with nanosized magnetic particles em- 
bedded in a polymeric matrix. The PANI phase surrounding cit-MAG nanoparticles prevents particle clustering 
leading to a structure with some regularity in terms of particle shape and position. In this binary material system, 
the sample can be idealized as an array of isolated spheres, with a well-defined size distribution, packed in a local 
HCP order and presenting long-range disorder. The pattern observed in the film cross-sectional micrograph resembles 
those in samples consisting of microphase separated diblock copolymers, in which one of the polymer subchains is 
much longer than the other one. This similarity motivated us to select the CDS model for the phase segregation of 
block copolymers to simulate the morphology of the as-produced nanostructured magnetic films. The CDS technique 
was originally proposed by Oono and Pun^'^ and has been successfully employed to simulate the spatial structure of 
block copolymers^ with great computational efficiency. Let us remind that the polymeric matrix itself surrounding 
the maghemite nanoparticles we study in this paper has no relation with the diblock copolymers that served as an 
inpiration for the use of the CDS technique. It is only the final pattern generated by the CDS numerical algorithm 
that resembles the (PANI/cit-MAG)^ nanofilms, thus motivating us to use this specific technique. 

The CDS technique proposes models discrete in time and space. It is based on the sequential iteration of diffusely 
coupled maps that mimic the local segregation dynamics in terms of a local order parameter, ?/'(r, t), providing the 
number density difference between the two segregating species at the position f and time t. Non-segregated regions 
correspond to cells in which -0 ss 0. Segregated domains have either ip > or ip < 0, depending upon the predominant 
species. In the present case we have nanoparticles with fixed sizes, so we must seek a model in which the final 
segregated pattern has the same characteristic. Although the model provides the time evolution of the segregation 
process, we are interested only in the final pattern consisting of roughly spherical particles placed in a matrix with 
some local order. The iteration procedure converges to this pattern very quickly, so the method is extremely efficient. 

In the present study we consider a three dimensional lattice with periodic boundary conditions along all directions. 
The order parameter was defined such that he final pattern had spherical domains with ip < (maghemite nanopar- 
ticles) in a background with > (polymeric matrix). The system was initialized with ip{r,t) = "00 + Q!(r), where 
■00 is the off-criticality and a is a random number in the interval [-0.005,0.005]. The order parameter at time i + 1 is 
computed from the previous one according to Eq. [ij 

0(r, t + l) = 0(f, t) - V2 [/(f, t)] +B[i;,- ^(f, t)] , (1) 

where 



I(f,t) = Atanh[i/;(r,t)] - %p{r,t)+T^'^ ['^{r,t)] 



(2) 
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Here, A > 1 is the segregation strength and defines the maximum absolute value of ip. F is a phenomenological 
diffusion constant and the parameter i? <C 1 is responsible for the partial segregation. The symbol stands for the 
three-dimensional isotropic discrete Laplacian defined ad^ 
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where the subscripts nn, nnn and nnnn stand for "nearest-neighbor" , "next-nearest-neighbor" and "next-next-nearest- 
neighbor" , respectively. 

The spatial structure to be analyzed was simulated by iterating 2000 times Eqs. [l] and [2] on a 300 x 300 x 30 slab, 
with the 300 x 300 planes parallel to the xy plane. In all simulations we used the parameter values that are know to 
generate the pattern we are looking io^^. ipo — 0.3, A — 1.3, F — 0.5, and B = 0.02. 



B. Analysis of the simulated morphology 

After the CDS simulation one has the values of order parameter for all cells. The next step is to convert this 
information into Np sets of parameters {fi,rhi,Kf ,ei} defined as position (r^), total magnetic moment {frii), shape 
anisotropy constant (Kf) and easy magnetization axis (ci) for each nanoparticle. 

The first step is to identify the particles. For that, a three-dimensional extension of the Hoshen & Kopelman 
algorithm^ for cluster identification and labeling was used. This analysis lead to the total number of clusters {Np) 
and to the number of cells belonging to each cluster {{ui}). Afterwards, we determined the clusters positions by 
calculating the center of mass of the group of cells forming each cluster. Next comes the analysis of the particles shape 
for determination of shape anisotropy properties (easy axis direction and strength) of each cluster. We considered 
that the particles had a prolate ellipsoidal shape. In this case the shape anisotropy constant could be calculated as 
foUowpl 
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where a and c are lengths of the shorter and longer particle axis, respectively. The ratio r = c/a and the easy axis 
direction can be determined by diagonalizing the inertia matrix for each particle. 

Another issue is the calibration of the length scale. The pattern generated by the CDS model has an arbitrary 
length scale. In order to obtain a good quantitative comparison with experiment we calibrated the simulated pattern 
by defining that the simulated average diameter (i^cDs) had the same value as in the real sample (I^tem); that is 
DcDS — 7.5 nm. With this, the length of each cell is 



'cell 



TEM 



(6) 



where ric is the average number of cells in the clusters. The magnetic moment of each cluster depends on its 
number of cells (n;) and can be obtained from the experimental value of the saturation magnetization (Ms) through 



For this simulated system, the average shape anisotropy constant was (if*) = 7.2 x 10 erg/cm^ as 



calculated through Eq. [ij The simulated particle size distribution in the ensemble is presented in Fig. ^h). 



C. Monte Carlo simulation 



In order to proceed with the simulation of the magnetization curves it is necesary to specify which energy terms will 
be taken into account in the MC simulation. According to our simulated ensemble, the distance i?nn between grains 
was relatively large, so that inter-particle interaction, such as direct exchange and RKKY, were neglected. Within 
this approximation, the magnetic dipolar interaction was the only considered one among particles. We wrote down 
this interaction in terms of magnetic dipoles positioned at the center of the particles. Finally, the anisotropy energy 
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included shape and crystalline contributions. Including the terms mentioned above we defined the total energy of the 
Np particles in the presence of an external magnetic field [H) as: 



1=1 



E 

j>i 



rhi ■ rhj - 3(TOi • hj)irnj ■ fij) 



(7) 



where the first term denotes the particles Zeeman term and accounts for shape and crystalline anisotropy energies. 
Based on the experimental data available from the literature^^lSl] -^^e have assumed uniaxial crystalline anisotropy for 
cubic maghemite, thus using the value Ku = 1.25 x 10^ erg/cm^ in the simulations. Note that the third term in Eq. [t] 
refers to the particle-particle dipolar interaction, in which fij defines the distance between particles i and j. In order 
to mimic the nanocompositc film we have applied periodic boundary conditions for the magnetostatic interaction on 
the xy plane. 

The MC simulations were performed using the Metropolis algorithm. In short, a new orientation for the magnetic 
moment of a randomly chosen nanoparticle is set, then the energy difference (A_E) between old and new orientations is 
calculated and accepted with probability p = min [1, exp(— AiS/fceT)]. The new orientation for the magnetic moments 
is chosen according to an algorithm analogous to the one proposed by Nowak et al^, in which a randomly oriented 
vector V with modulus G|TOi| is added to the moment m^. By changing G one can control how different new and old 
directions might be. With this algorithm the magnetic moment is restricted to a cone centered in with an opening 
angle (?niax = sin~^(G) (see Fig. jsj. We have used ^max = 15° throughout the simulation. One MC step (MCS) 
corresponds to examining, in average, all the magnetic moments. 



IV. RESULTS AND DISCUSSION 

Figure [l] shows the experimental (a) and simulated spatial structures (b). A qualitative analysis shows good 
agreement between both structures. Following with the morphology comparison, we calculated the size histogram 
that was well fitted by a log-normal distribution with (t = 0.11 nni. Simulated and experimental histograms are shown 
in Fig. [2] and the comparison between them reinforces the adequacy of the simulated sample. Finally we calculated 
the pair correlation function, g{r), from which one can estimate the average distance between nearest-neighbors (i?nn)- 
Figure [4]jb) displays g{r) and clearly indicates short-range order at distances smaller than 30 nm. From the first g{r) 
peak position we concluded that the average center-to-center nearest neighbor distance is i?nn = 10 nm. Considering 
that the average particle diameter is 7.5 nm, the average interparticle separation is about 2.5 nm. In the real sample 
this number represents the average distance of neighboring particles, that is a surface-to-surface spacing of about 2.5 
nm, which is filled by the hosting polymeric material, in very good agreement with the observed TEM cross-sectional 
micrographs in Fig. |4](a). 

The simulated ZFC/FC curves were obtained starting from a demagnetized state at 2 K. The temperature was then 
increased at a rate of = 50 MCS/2K, up to 300 K, under an applied in-plane magnetic field of 36 Oe. Here MCS 
stands for Monte Carlo steps. For the simulation of the FC curve we started the simulation with the last configuration 
of the ZFC curve and the temperature was decreased down to 2 K at the same rate as in the ZFC curve. This 
procedure was repeated 2000 times and the results were averaged out over all statistically equivalent realizations. For 
simplification we have considered = (K^) (its average over the ensemble). Experimental and simulated ZFC/FC 
curves are presented in Figures [sf^ a) and[5jb), respectively. The values of Tmax in the ZFC curves were 54 and 53 
K for the experimental and simulated curves, respectively. Since Tmax scales with the blocking temperature Tb we 
can draw conclusions about the amount of interaction in the samples through the variation of Tmax- The decreasing 
of £ (the number of bilayers) in the (PANI/cit-MAG)^ magnetic nanocompositc films, causes the reduction of Tmax 
from 54 K (£ = 50, 25) to 45 K (£ = 10). This is consistent with the picture of particle interaction effects reported in 
the literatur^^ as the content of the magnetic nanoparticle incorporated within (PANI/cit-MAG)^ reduces with the 
reduction of £. The small discrepany in Tmax from 53 K (simulated) in comparison to 54 K (experimental) , as shown 
in Figures |5j indicates that the two samples (simulated and experimental) are indeed quite similar, though the real 
sample might present stronger particle-particle interaction than the simulated one. Interesting to note that Tmax = 54 
K was observed in both nanocompositc films {i — 50, 25), indicating that the nanoparticle packing has reached a limit 
and the only difference between the two samples is actually the film thickness. Furthermore, while comparing with 
the simulated ZFC/FC curves the observed fiattening of the experimental ZFC/FC curves {£ = 10,25,50) indicates 
the mean-field effect, more likely due to the increasing of the samples demagnetizing factor as a consequence of the 
particle chain formatiorP^l. While the simulated ZFC/FC curves shown in Fig. ^h) were derived from a system 
comprising an ensemble of single nanoparticles, though including dipolar interaction among them, as described in Eq. 
[7] the experimental data do indicate that particle clustering takes place in the real samples we investigated here. 
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The discussion regarding the data provided by the ZFC/FC curves is corroborated by the experimental (£ = 25) 
and simulated hysteresis curves shown in Figures [6]j a) and|6][b), respectively. Simulations of hysteresis cycles were 
performed at T = 2 and 300 K while applying an in-plane magnetic field. For the simulated cycles the saturation field 
Hs was determined as the field at which the magnetization reaches 99% of the sample's saturation magnetization. After 
reaching saturation the in-plane applied field was swept between +Hs and — i?s at a rate oiv = f 00 MCS/10 Oe. Notice 
the similarity of the room temperature (T = 300 K) reduced remanence values for the experimental {Mr /Ms — 0.54) 
and simulated [M/Mg = 0.49) curves, shown in Figure |6] Additionally, notice the similarity of the experimental and 
simulated hysteresis cycles at T = 300 K, with minimal remanence (M^j) and coercivity [He)- Nevertheless, the 
low temperature (T = 2 K) reduced magnetization, remanence and coercivity values observed in the data provided 
by Figures [6j^a) and [6]jb) are quite different. We found the experimental and simulated normalized magnetization 
{M* /Ms) at H* = 1 kOe equal to 0.61 and 0.93, respectively. This is a clear indication that particle clustering takes 
place in the real (PANI/cit-MAG)25 sample. In support to this model picture we found the experimental normalized 
magnetization {M* /Ms) at H* = 1 kOe equals to 0.68 for the (PANI/cit-MAG)io sample, indicating reduced particle 
clustering as the particle concentration reduces (from I ~2b down to ^ = 10) within the magnetic nanocomposite film. 
Additionally, from Figure [6] we found the experimental and simulated low temperature (T = 2 K) reduced remanence 
{Mji/Ms) equal to 0.32 and 0.53, respectively. The simulated value of Mji/Ms = 0.53 is very much closer to the 
value predicted by the Stoner-Wohlfarth theory {M^/Ms — 0.50) of non-interacting uniaxial single domain particles 
at r = kE21. 

On the other hand, particle-particle interaction is expected to reduce the values of M/f/M^^^, as observed in the 
nanocomposite films investigated here. In agreement with the discussion above we found for the (PANI/cit-MAG)io 
sample Mr/Ms = 0.35 at T = 2 K, confirming the reduced particle content within the magnetic nanocomposite 
film of this sample in comparison with the (PANI/cit-MAG)25 sample. Likewise, the analysis of the coercivity (He) 
values obtained from the data shown in Figure |6] reveals the expected tendency of enhancement with reduction of the 
particle-particle interactiorP^. The low temperature (T = 2 K) values of coercivity observed in the experimental and 
simulated curves shown in Figures [6]ja) and[6f^b) were He = 210 and He = 274 Oe, respectively. 

A final comment regarding the success of the simulation approach developed herein is based on the correlation 
function presented in Fig. [4] From Fig. [4]d wc found the system is order up to about 30 nm in length, indicating that 
we should expect a tendency of saturation of the magnetic properties which are dependent upon the particle-particle 
interaction in nanocomposite films (like ours) with thicknesses higher than about 30 nm. Recently, we found the 
thickness of the (PANI/cit-MAG)^ system linearly increasing with £, the slope increasing with the concentration of 
the suspended nanoparticles within the magnetic fluid sample used for fabricating the film^^S] Indeed, we observed 
Tmax, Mb,/ Ms and He systematically changing as one goes from sample (PANI/cit-MAG)io to sample (PANI/cit- 
MAG)25. Nevertheless, only slight changes in all these values (Tmax, M^JMs and He) were observed as one goes from 
sample (PANI/cit-MAG)25 to sample (PANI/cit-MAG)5o, indicating that above £ = 25 the film reached a thickness 
higher than the correlation length of about 30 nm. These findings represent a strong support to the usefulness of the 
simulation approach presented here. 

V. CONCLUSION 

In this study nanofilms of citrate-coated maghemite nanoparticles embedded in polyaniline, (PANI/cit-MAG)^, 
were produced by the layer-by-layer technique, experimentally investigated via high resolution microscopy (HRM) 
and magnetization measurements (ZFC/FC and M x H curves) and employed as a model system to evaluate the 
effectiveness of the simulation approach herein proposed. We found the HRM data (samples morphology) in very 
good agreement with the results provided by the Cell Dynamic System (CDS) simulation. For a typical sample the 
simulated distribution function, g{r), provided the first peak position at 10 nm for nanosized particles with 7.5 nm 
in diameter, thus leaving about 2.5 nm for the surface-to-surface distance. This simulated value is in very good 
agreement with the surface-to-surface distance of 3 nm obtained from the HRM micrograph. 

The samples morphology provided by the CDS simulation was then used as input to Monte Carlo (MC) simula- 
tions of ZFC/FC curves and M x H loops, including energy terms corresponding to shape and crystalline magnetic 
anisotropy and dipolar interactions. The simulated curves presented good correspondence in several important fea- 
tures regarding real samples, providing good estimation and understanding for the ZFC curve maximum (Tmax), 
remanence (MfjjMs), and coercivity (He)- For a typical sample we found the simulated Tmax (53 K) in excellent 
agreement with the experimental data (54 K). Nevertheless, for a typical sample we found experimental values of 
remanence and coercivity sistematicaly lower {Mr/Ms = 0.32 and He = 210 Oe) than the values we found in the MC 
simulation {M^/Ms = 0.53 and He = 274 Oe). We understood that this difference is due to the particle clustering in 
the real samples, which is expected to lower both M^/Ms and He with respect to the values for a system consisting 
of isolated nanoparticles. 
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In summary we conclude that the combination of CDS and MC simulations provided a realistic description for the 
as-produced nanofilms. The method developed herein can be employed to design optimized magnetic nanocomposites 
and to predict their respective properties with great accuracy. This novel approach would also find application in 
fields beyond that investigated here, spanning from biology and medicine to engineering and nanoelectronics, where 
phase segregated systems play a major role. 
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FIG. 1: Typical surface morphology of a (PANI/cit-MAG)io film (a) and simulated (b) cross-section structure after CDS 
modeling. 
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FIG. 2: (a) Particle size histogram of the cit-MAG colloidal sample, obtained from the TEM images. The solid line results 
from the curve-fitting of the data using the log-normal distribution function, (b) Particle size histogram of the average particle 
diameter after CDS simulation. 
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FIG. 3: At each trial step a magnetic moment rhi (black) is selected. A random vector v with modulus G\fhi\ (red) is added to 
the latter. The resulting vector v' (blue) gives the direction of the new magnetic moment. After rescaling it in order to keep 
the same modulus as rhi, it will be the new magnetic moment m'i We can control the angle 6 between fhi and m'i by choosing 
the appropriate G. 




FIG. 4: Typical cross-section TEM micrograph of a (PANI/cit-MAG)io film (a); and simulated pair distribution function, g{r), 
for the 300 x 300 x 30 grid (b). 
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FIG. 5: Experimental (PANI/cit-MAG)25 nanofilm (a) and simulated (b) ZFC/FC curves obtained at a steady field of 36 Oe. 
Solid lines on (b) are guide to the eyes only. 




FIG. 6: Experimental (PANI/cit-MAG)25 nanofilm (a) and simulated (b) hysteresis curves at different temperatures. Solid 
lines on (b) are guide to the eyes only. 



